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^ I Abstract: The understanding of the origin of dark matter has great importance for cosmology 

^ I and particle physics. Several interesting extensions of the standard model deaUng with solution 

Q I of this problem motivate the concept of hidden sectors consisting of SU{3)c x SU{2)l x U{1)y 

00 ■ singlet fields. Among these models, the mirror matter model is certainly one of the most interest- 

. ■ ing. The model explains the origin of parity violation in weak interactions, it could also explain 

Q ■ the baryon asymmetry of the Universe and provide a natural ground for the explanation of dark 

O . matter. The mirror matter could have a portal to our world through photon-mirror photon mixing 

(e). This mixing would lead to orthopositronium (o — Ps) to mirror orthopositronium oscillations, 
the experimental signature of which is the apparently invisible decay of o — Ps. In this paper, we 
/S ' describe an experiment to search for the decay o — Ps—?' invisible in vacuum by using a pulsed slow 

C^ ■ positron beam and a massive 4k BGO crystal calorimeter. The developed high efficiency positron 

tagging system, the low calorimeter energy threshold and high hermiticity allow the expected sen- 
sitivity in mixing strength to be £ ~ 10^^, which is more than one order of magnitude below the 
current Big Bang Nucleosynthesis limit and in a region of parameter space of great theoretical and 
phenomenological interest. The vacuum experiment with such sensitivity is particularly timely in 
light of the recent DAM A/LIBRA observations of the annual modulation signal consistent with a 
mirror type dark matter interpretation. 

Keywords: Positronium; Dark matter; Hidden sectors. 



Contents 



|l[ Introduction 



Mirror Matter model 

|2.l| Higgs portal into mirror world at the CERN LHC 

[2.2| Positronium portal into mirror world 

[2.3| Experimental evidence for Dark Matter of Mirror type 



i 


Experimental technique 


i 


The setup 




1 1 The slow positron beam 




12 o — Ps production target 




4.3 The vacuum cavity 




4.4 Positron tagging system 




|4.4. 1 Trigger efficiency and confidence level 




[4.5 The photon detector 



Background estimation 

IsTT] Fast backscattered o 



Ps from the carbon foil 



5^ Backscattered positrons 

Positron backscattering from the carbon foil 
Positron backscattering from the SiOj target 



5.2.1 



5.2.2 



^ Sensitivity 
^ Summary 





m 

13 
14 
B 

23 



23 
25 
26 
29 
30 



35 



1. Introduction 

Cosmological observations of galactic rotational curves ^ and the gravitational lensing [Q, ^] give 
strong evidence for the existence of dark matter (see e.g. [Q] for an excellent review). In addition, 
the recent confirmation by the DAMA/LIBRA experiment ^ of the annual modulation signal ob- 
served by the DAMA/Nal [|6|] seems provide the first direct terrestrial experimental observation of 
the existence of non-baryonic dark matter in our galactic halo. The need to explain these hints of the 
presence of dark matter provides one of the strongest indications for the existence of new physics 
beyond the Standard Model (SM) and the identification of the origin of dark matter is a task of 
enormous importance for both particle physics and cosmology. At present, the most popular candi- 
dates for the (thermal-produced) dark matter are the so-called weakly interacting massive particles 



(WIMPs), which are e.g. lightest supersymmetric particles, Kaluza-Klein particles in universal ex- 
tra dimension models, axions etc... However, despite of significant efforts the experiments looking 
for WIMPs lead so far to negative results, thus, pushing further possible WIMP searches into a very 
high-energy and/or high sensitivity frontiers, for a recent review see e.g. [^ and references therein. 

An additional natural ground for the explanation of dark matter is provided by a class of in- 
teresting theoretical models introducing the concept of "hidden" sectors consisting of SU{3)c x 
SU{2)l X U{1)y singlet fields. It is worthwhile to note that, even in the SM, some of the matter 
fields are singlet under one or more of the colour and electroweak gauge groups. Thus, the exten- 
sion to include a further sector which transforms under the new but not under the familiar gauge 
symmetries is not particularly exotic from a theoretical viewpoint. 

The sensitivity of searches for the new singlet particles depends in detail on their couplings and 
mass scale, e.g. if the mass scale of a hidden sector is too high, it is experimentally unobservable 
and indeed is hidden. Then the question arises: could the important sensitive searches for the 
hidden sectors be performed at low energy frontier ? The answer for this question is definitely 
positive. For example, there is a class of models with at least one additional U{1) gauge factor 
where the corresponding hidden gauge boson could be light, or even massless ||8|, ^, 10 1, for a 



recent review see [11]. With respect to the mass and lifetime range for hidden sector states, the 
niche with masses in the range of a few MeV to a few GeV, with lifetimes of less than 1 second, 
is also quite intriguing as many of the most severe astrophysical and cosmological constraints that 
apply to lighter states are weakened or eliminated, while those from high-energy colliders are often 



inapplicable [12] 



Recent hints from terrestrial and astrophysical anomalies, in particular first CDMSII result 
[13] and annual modulations from DAMA/LIBRA at Gran Sasso [^, ^], positron and/or electron- 
positron pair excess in cosmic ray from PAMELA [ |l4| ] and ATIC [ |l5| ] (however, the ATIC excess 
seems has not been confirmed by Fermi -LAT [|^]) are found to be consistent with an interpretation 
in context of dark matter charged under new gauge fields [[T7|]. Thus, if those observations are 
really signals for DM they would indicate particularly exotic forms of DM. This suggests that 
testing hidden sectors with the relatively light mass scale states in the MeV-to-GeV range weakly 
coupled to the SM model represents an intriguing possibility. 

Then it is important to ask which hidden sector models could be tested with the best experi- 
mental sensitivity? Certain classes of hidden sectors may be tested, but it is clear that not many 
generic scenarios, which are viable, could be strongly constrained and even excluded on the base 
of precision measurements. Among numerous candidates that have been discussed, one of the 
most promising, which could reconcile the DAMA annual modulation signal with the results of the 



CDMS II experiment, is mirror type dark matter [llq, 19]. In addition, the CoGeNT collaboration 



has recently reported a rising low energy spectrum in their ultra low noise germanium detector 



[20]. This is particularly interesting as the energy range probed by CoGeNT overlaps with the 



energy region in which DAMA has observed their annual modulation signal. It has been recently 



shown [21], that the mirror dark matter candidate can simultaneously explain both the DAMA an- 
nual modulation signal and the rising low energy spectrum observed by CoGeNT. This constitutes 
a model dependent confirmation of the DAMA signal and adds weight to the mirror dark matter 
paradigm. 

In this paper, we show that the proposed search could result either in a strong evidence of the 
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existence of dark matter of the mirror type, or completely exclude this type of new physics beyond 
the standard model. It should be noted that mirror baryons are naturally dark, stable and massive. 
Currently, it seems that this type of matter could also explain in a natural way the coincidence 
between visible and dark matter densities in the universe (Hg = 0.044 and ^dm = 0.26) [ p2| , ^ ]. 
Interestingly, it could also provide an explanation of the controversial behaviour in galaxy cluster 



collisions [24|. Furthermore, if mirror matter is present in our universe it would mean that Parity 
(spatial-inversion) is an unbroken symmetry of nature. This gives an exciting motivation for testing 
this model by a laboratory experiment at low energies. 



The rest of the paper is organized as follows. In Sec. |2.3| , we review the mirror model and 
work out briefly the mirror Higgs boson physics at the CERN LHC, and in more details, the effect 
of oscillation of ordinary positronium to the mirror one and its experimental consequences. We also 
present the calculations we performed in order to estimate the effect of matter and external fields 
on the oscillation probability. In Section |3[ we present the experimental technique. In Section 0, 
the design of the experiment and the detector components are described in detail. The simulations 
of the signal and background sources, as well as the expected sensitivity are discussed in Section || 
and ^, respectively. Section ^ contains concluding remarks. 

2. Mirror Matter model 



Mirror matter was originally discussed by Lee and Yang [ ]25| ] in 1956, after their discovery of 
parity violation (for an excellent recent review on this subject see [26]). In order to save parity 
conservation they suggested that the transformation in the particle space corresponding to the space 
inversion x — )• — x should not be the usual transformation P but PR, where R corresponds to the 
transformation of a particle (proton [^]) into a reflected state in the mirror particle space. After 
the observation of parity non-conservation. Landau assumed [^] that R=C, i.e. he suggested to 
identify antiparticles with the mirror matter but then CP must be conserved, which we know is not 



the case. The idea was further developed by A. Salam |28], and was clearly formulated in 1966 as 
a concept of the mirror universe by Kobzarev, Okun and Pomeranchuk [^]. In their paper, they 
have shown that ordinary and mirror matter can communicate predominantly through gravity and 
proposed that the mirror matter objects can be present in our universe. 

Since that time, the concept of mirror matter has found many interesting applications and 
developments. In the 80's, it has been boosted by superstring theories with Eg, x E^ symmetry, 
where the particles and the symmetry of interactions in each of the Eg, groups are identical. Hence, 
the idea of mirror matter can be naturally combined in these models [^]. 

Nowadays, mirror matter models exist in two basic versions. The symmetric version, proposed 
earlier, was further developed and put into a modern context by Foot, Lew and Volkas [|1|]. The 
asymmetric version was proposed by Berezhiani and Mohapatra [^]. 

In the following we will concentrate on the symmetric model since it is the most interesting 
from a dark matter perspective and it could provide, as we will see, an experimental signature 
related to positronium. In the symmetric mirror model, the idea is that for each ordinary parti- 
cle, such as the photon, electron, proton and neutron, there is a corresponding mirror particle of 
exactly the same mass as the ordinary particle. R-parity interchanges the ordinary particles with 
the mirror particles so that the properties of the mirror particles completely mirror those of the 
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ordinary particles. For example, the mirror proton and mirror electron are stable and interact with 
the mirror photon in the same way in which the ordinary proton and electron interact with the or- 
dinary photons. The mirror particles are unlikely to be produced in laboratory experiments just 
because they couple very weakly with the ordinary particles. In the modem language of gauge the- 
ories, the mirror particles are all singlets under the standard G = SU{3) ® SU{2)i U{1)y gauge 
interactions pT|]. The mirror particles interact with a set of mirror gauge particles, so that the 
gauge symmetry of the theory is doubled, i.e. the minimal gauge group of the new mirror model 
is GsM <E) G'sM = SUi3)c «> SU{2)l (S)U{1)y'S) SU'{3)c «> SU{2y^ (g) ^'(l)y (the ordinary particles 
are, of course, singlets under the mirror gauge symmetry) [jsTl]. The gauge quantum numbers under 
GsM ^ G'^M foJ" ^^^ usual and new fermion fields are 

(Li)'~(l,2,-l)(l,l,0), (4y~(l,l,0)(l,2,-l) 
(^«)'~(1,1,-2)(1,1,0), (eiy~(l,l,0)(l,l,-2) 

(QlY ~ (3,2, i)(l, 1,0), {Q'^y ~ (1, 1,0)(3,2, i) 

(urY ~ (3, 1, ^)(1, 1,0), {u'J ~ (1, 1,0)(3, 1, ^) 

(Jfi)'~(3,l,-^)(1,1,0), (4)'~(1,1,0)(3,1,-^) 

(2.1) 

with i the family index. In the left-right symmetric models parity is extended to a new type Z2 
discrete symmetry which transforms the left-handed field to the right-handed one for the same 
fermion. However, the new Z2 parity symmetry that we can define now, is 

x^-x,t^t,G''^ G'^, W ^W^, 
B^ o B'^, U o 4, eR o el, Ql o Qr, 

ur o m^, cIr -f-)- d'l^ (2.2) 



We see, that under Z2 of ( |2.2| ), the left-handed sector of the fermion field can transform to the right- 
handed sector of a different fermion field, namely, the mirror fennion field. Thus, the Z2 symmetry 
can be interpreted as a parity symmetry {x — ;■ —x), if the roles of left and right chiral fermion fields 
are interchanged in the mirror sector. Parity is conserved because the mirror particles experience 

V +A (i.e. right-handed) mirror weak interactions while the ordinary particles experience the usual 

V —A (i.e. left-handed) weak interactions. 

An exact Lagrangian Z2 symmetry interchanging ordinary and mirror particles is hypothesized, 
which means that all the couplings in the mirror sector are the same as in the ordinary sector. While, 
the ordinary and mirror particle sectors can interact between each other in a number of ways. The 
first is through gravitation, with immediate consequences for the dark matter problem and astro- 
physics [p3|]. Non-gravitational interactions can be induced through the mixing of colourless and 



neutral particles with their mirror counterparts. Neutrinos [|34|], the photon [ [35[ , |36| , |370 , and the 



physical neutral Higgs boson [ |38[ , |39p can mix with the corresponding mirror states. Coloured 
and/or electrically charged particles are prevented from mixing with their mirror analogues by 
colour and electric charge conservation laws. 



It is known pl[ ] that there are two renormalizable and gauge invariant Lagrangian terms cou- 
pling the ordinary and mirror sector together: U {1)y — U {1)'y gauge boson kinetic mixing, and 
ordinary-mirror Higgs scalar interactions, i.e. 

f^^v^'" (2-3) 

and 

^vImVsmvIvh (2.4) 

Based on this observation, one can extend the SM by doubling the ordinary fermion, gauge, 
and Higgs fields [|1|]. Thus, the new mirror fermions are natural singlets of the SM gauge group, 
and they (nucleus if there exists mirror SU{3)c) can be the candidates for dark matter. 

Mirror matter is invisible to us because it does not interact with ordinary photons and naturally 
constitutes a dark matter candidate. One should stress that the fact that our and mirror-sectors 
have the same micro-physics, does not imply that their cosmological evolutions should be the 
same. Indeed, if mirror particles had the same temperature in the early universe as ordinary ones, 
this would be conflict with Big Bang Nucleosynthesis (BBN). The BBN limit on the effective 
number of extra neutrinos implies that the temperature of the mirror sector T' must be at least 
about twice smaller than the temperature T of the ordinary sector allowing mirror baryons to be 
a viable candidate for dark matter. In particular, the mirror dark matter scenario would give the 
same pattern of Cosmic Microwave Background and Large Scale Structure as the standard CDM if 



T' /T < 0.2 [^2, 33, 40]. In addition, the baryon asymmetry of the Universe can be generated via 



out-of- equilibrium B — L and CP violating processes between ordinary and mirror particles [41] 



whose mechanism could explain the intriguing puzzle of the correspondence between the visible 



and dark matter fractions in the Universe, naturally predicting the ratio Q-um/^b — 5 []42|]. 

2.1 Higgs portal into mirror world at the CERN LHC 

In the mirror model every standard particle, including the physical neutral Higgs boson, is paired 



with a parity partner. The interaction of Eq. ( |2.4| ) and unbroken parity symmetry forces the mass 



eigenstate Higgs bosons to be maximal mixtures of the ordinary and mirror Higgs bosons: 

H^ = ^ (2.5) 

V2 

Each of these mass eigenstates will therefore decay 50% of the time into invisible mirror particles. 
The total decay rate of //+ or H is the same as that for a SM physical neutral Higgs boson of the 
same mass. This may result in dramatic consequences for the LHC, making the significance of the 
Higgs signal at LHC lower due to decreasing of the Signal/Background ratio if the mass splitting 
is large compared to the Higgs mass resolution at the LHC. Note also that each mass eigenstate 
couples to ordinary particles with strength reduced by 1/2 compared to the coupling of the standard 



Higgs boson to those same particles [38] 



Double Higgs peak observation for Eq. ([2.5D, would give a clear and interesting signature 



for the Large Hadron Collider(LHC) which could thus establish the existence of the mirror world. 
However, for this effect to be observable the mass difference between the two eigenstates must 
be sufficiently large. The cosmological constraints from Big Bang Nucleosynthesis on the mass 



difference parameter have been studied in [38]. To summarize, one can see that the Higgs sector 



may play an important role in detecting mirror particles, which can be the candidates of dark matter 



and appear as missing energy in the detectors at the LHC [ )38| , [39| ]. Another intriguing possibility 
discussed in the next Sec. is related to the possible observation of orthopositronium to mirror 
orthopositronium oscillations. 

2.2 Positronium portal into mirror world 

Positronium (Ps), the positron-electron bound state, is the lightest known atom, which is bounded 
and self-annihilates through the same, electromagnetic interaction. At the current level of experi- 



mental and theoretical precision this is the only interaction present in this system [|430. This feature 
has made positronium an ideal system for testing the accuracy of QED calculations for bound states, 
in particular for the triplet (PS'i) state of Ps, orthopositronium (o — Ps). Due to the odd-parity un- 
der C-transformation, o — Ps decays predominantly into three photons. Due to the phase-space and 
additional a suppression factors, as compared with the singlet (1^5o) state (para-positi^onium), the 



'slowness" of o — P^ decay rate gives an enhancement factor ~ 10 in sensitivity to an admixture 



of new interactions which are not accommodated in the Standard Model [44, 45]. Glashow realized 
that the orthopositronium system provides one sensitive way to search for the mirror matter [p^]. 
Glashow 's idea is that if a small kinetic mixing between ordinary and mirror photons exists [^], it 
would mix ordinary and mirror orthopositronium, leading to maximal orthopositronium - mirror 
orthopositronium oscillations (see Fig. |l]). Since mirror o-Ps' decays predominantly into three mir- 
ror photons these oscillations would result in o — Ps — )• invisible decays in vacuum. Photon-mirror 
photon kinetic mixing is described by the interaction Lagrangian density 

L = eF^''F'^^, (2.6) 

where F^^ {PLv) is the field strength tensor for electromagnetism (mirror electromagnetism). 

Together with the Higgs- mirror Higgs quartic couple X^^ ^'^' , these are the only renormal- 
izable and gauge invariant terms that can be added to the SM Lagrangian. The effect of ordinary 
photon - mirror photon kinetic mixing is to give the mirror charged particles a small electric charge 
[ ^Ij ^ ^5|]. That is, they couple to ordinary photons with chai^ge 2ee^ . 

Orthopositronium is connected via a one-photon annihilation diagram to its mirror version (o- 
Ps') [^]. This breaks the degeneracy between o-Ps and o-Ps' so that the vacuum energy eigenstates 
are {o — Ps + o — P^)/^/2 and {o — Ps — o — Ps')/V2, which are split in energy by 

AE = 2hef, (2.7) 

where / = 8.7 x 10^ MHz is the contribution to the ortho-para splitting from the one-photon 
annihilation diagram involving o-Ps [^]. Assuming a mixing strength of 4 x 10^^ (as suggested by 
the DAMA results), one obtains an energy splitting of AE = 2.9 x 10^'^ eV. Thus, the interaction 
eigenstates are maximal combinations of mass eigenstates which implies that o-Ps oscillates into 
o-Ps' with a probability: 

P{p-Ps^o-Ps')=sm^Q.t, (2.8) 



'Note, that the direct experimental bound on e from searches for 'milh-charged' particles is e ~ 10 
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Figure 1. The double degeneracy between orthopositronium mass eigenstates of ordinary (o-Ps) and mirror 
(o-Ps ') is broken when a small mixing term is included. 



where D. = 2nef. 

The simplest case ofo — Ps^o — Ps' oscillations in vacuum [ 35 ] leads to an apparent increase 
in the decay rate because the mirror decays are not detected. The number of o-Ps N satisfies 

: cos^ Q.t ■ e^^'"' ~ exp[-t{TsM + ^^0] > 



A^: 



(2.9) 



where Tsm is the Standard Model decay rate of o - Pi [|7[ |8|, |9|]. Thus T'^ff ; 
leads to a branching ratio of: 



:r5M(i+n7r 



SMI 



Br{o — Ps^ invisible) 



2{2Keff 



(2.10) 



Yl^ + A{2Kefr 

The above calculation is not applicable to an experiment performed with a cavity confining 
the positronium, because in this case the collision rate is not zero and the loss of coherence due to 



the collisions must be included in the calculation [36, 37]. In those papers, the collision rate was 



assumed to be much larger than the decay rate Tcoii ^ ^sm [37], thus, an approximate solution for 
this case was found. In the proposed experiment, we are interested to see the effect on the oscillation 
probability for 1-2 collisions per lifetime. Furthermore, external electric or magnetic field induce 
additional splitting of o-Ps and o-Ps' states affecting the probability P{o — Ps^o — Ps'). 

Only the second order Stark shift contributes to the positronium in the ground state. The 
energy shift can be calculated in the same way as hydrogen and is given by [BO] : 



^stark 



1 „2 



(2.11) 



where 4.54ao = ^<^q is the polarizability of o-Ps (ao = 0.1 nm is the Bohr radius of positronium), 
Co is the vacuum permittivity and E is the electric field. Magnetic fields affects only the triplet state 
with quantum number m=0. For the m=ib 1 states, the Zeeman effect is zero because the magnetic 
moments for positron and electron are opposite. The energy contribution for m=0 can be calculated 
with: 

t.2 ' 



Ar 



■ + 



Ah' 



)^ + hh\Bl 



(2.12) 



where A = 4.92 x 10^^ eV 's ^ is the hyperfine splitting constant of the ground state, g\ = 2.81 x 
jqIO x^^s^' the gyromagnetic ratio, Bq the magnetic field. 



To include these effects in the calculation, we use the density matrix approach. The Heisen- 
berg's equation of motion 

P = UH,p] + ^Prel (2.13) 

in at 

can be solved to find the density matrix 

pi^)=(''i\''i\] (2-14) 

that describes the evolution of the states with time. The hamiltonian H of the system is given by: 

On 
2 



"= -QS -«5l <''« 



where A = Astark+ ^zeeman is the shift of the ground state introduced by magnetic and electric fields. 
The effect of those fields on o-Ps' can be neglected because of the very weak coupling. The term 
0)12 is the splitting between the two vacuum eigenstates. To describe all the processes (decays 
and collisions) that return the ensemble in thermal equilibrium, i.e. destroy the coherence of the 
oscillation, we use the relaxation term p^/. Its derivative is equal to: 

d _( -riPnit) -(]l+ii + y„,„)pi2(0\ 

dtP-'-[-C-^ + r,^„)p,,{t) -np22{t) ) ^'-''^ 

where the constants /i and 72 are the decay rate in vacuum of oPs and oPs', respectively. The term 
Ycoii is the collision rate of o-Ps with the cavity walls. Substituting (2.15) and ( 2.16[ ) in Eq. (2.13), 
one obtains the four differential equations: 

pii(0 = -ripn{t)-U^npi2it)-^hp2i{t)) 

n 

Piiit) = -Y2p22{t)-U-Q.hpi2{t)+Q.np2i{t)) 
n 

Puit) = -{^^ + Ycou)Pl2{t) 

-^{^npn{t) + -conhpnit) + -{A + (Oi2)hpi2{t) - anp22{t)) 



P2l(0 = -(^ + 7co//)P2l(0 



-^{-^hpn{t)--0}i2fip2iit)--{A + 0)i2)fip2i{t)+anp22it)) (2.17) 



that can be solved numerically (we used Mathematica [ |51p for this purpose). At t=0 we start with 
a pure o-Ps state, thus, we set the initial conditions: Pii(O) = l,pi2(0) = 0,p2i(0) = 0,p22(0) = 0. 
We assumed that without the coupling the vacuum eigenstates of o-Ps and o-Ps' are degenerate and 
thus we set cou = in the calculation. To find the numerical solution, we also used 7i = 72 = 1/142 
ns^^ implying that o-Ps and o-Ps' have the same lifetime in vacuum. To find the branching ratios 
(BR) one has to integrate over t pn and P22, thus the BR is: 

JoPii(x)(ix + JoP22(x)(ix 



where the upper Umit of the integration can be chosen as the data acquisition gate for the calorime- 
ter. Assuming e = 4 x 10^^ as suggested by the DAMA/LIBRA results, one can calculate the 
oscillation probability (P22(0) ^^^ the branching ratios for o-Ps into o-Ps' for different values of 
the EB-fields. The results are shown in Figs. ^-^. As one can see, the branching ratio is not affected 
by the magnetic field (100 G) and by the electric field (less than 10 kV/cm) we are planning to use 
in the experiment. The effect of the collision is shown in ^ As expected, solving Eq. (2.18) for 



zero fields and no collisions leads to the same value of the BR obtained with Eq. (2. 10). 
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Figure 2. Oscillation probability and branching ratio for Br(o — Ps - 
netic field. The electric field was set to 5 kV/cm. 
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Figure 3. Oscillation probability and branching ratio for Br(o — Ps — ^ invisible^ as a function of the electric 
field. The magnetic field was set to 150 G. 



2.3 Experimental evidence for Dark Matter of Mirror type 

At present, there is some experimental evidence that mirror matter could exist, coming from cos- 



mology as well as from the neutrino physics [|52|]. Foot discussed implications of the DAM A 
experiment for mirror matter-type dark matter, which is coupled to ordinary matter through the in- 
teraction of Eq. (16) ^% p^ . It has been shown that the annual modulation signal measured by the 
DAMA/Nal experiment [^ can be explained by mirror matter-type dark matter if the photon-mirror 
photon mixing strength is in the region 



£ ~4x 10 



-9 



(2.19) 



Those results have been recently confirmed by the DAMA/LIBRA experiment with a 8.2a 
significance [9]. It was also pointed out 1 19] that the very recently published results of the CDM- 
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Figure 4. Branching ratio for Br(o — Ps — >■ invisiblejas a function of the number of collisions per lifetime of 
o-Ps with the cavity walls. The electric field was set to 5 kV/ctn and the magnetic field was 150 G. 
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Figure 5. The points are the annual modulation observed by DAMA/Nal. The line is the prediction of this 
modulation for mirror matter /|7q/. 



SII/Ge final exposure [13] with two events in the signal region with recoil energies close to the 
threshold can be explained by the interactions of mirror nuclei in the detector. Mirror matter is a 
promising candidate for dark matter since it can reconcile the null results from the other higher 
threshold experiments like CDMS/Si ||5|] and XENON 10 Q. Interestingly, this value of e is 
also consistent with all other known experimental and cosmological bounds, including SN1987a^ 



and the standard Big Bang Nucleosynthesis (BBN) bound []58[]. It was also confirmed that photon- 
mirror photon mixing of this magnitude is consistent as well with the more stringent constraints 
from cosmic microwave background measurements and large scale structure considerations [^]. 



It is also in the range of naturally small £-values motivated by grand unification models []22[]. 

If £ is as large as in Eq. (2.19|), the branching ratio Br{o — Ps ^- invisible) for the invisible 



decay of orthopositronium in vacuum can be found with Eq. ( [2.10 ) and is of the order: 

Br{o — Ps^ invisible) ~ 2 x 10^^. 



(2.20) 



The SN1987a limit e < 10 obtained in Ref. [pq] is actually much weaker. For a more detailed discussion of this 
and other constraints see Ref. |57^] 



-10- 



For comparison, the BBN limits [58] deduced from the successful prediction of the primordial ^He 
abundance are 

e < 3 X 10"^ (2.21) 

and 

Br{o -Ps-^ invisible) < 10"^ (2.22) 

respectively. 



Given the indications for the mirror world, coming from dark matter [53] and the neutrino 



physics anomalies [|52|, |60p, as well as the intuitive expectation that nature could be left-right sym- 
metric, it is obviously important to determine experimentally whether orthopositronium is a win- 
dow to the mirror world or not. 

3. Experimental technique 

The experimental signature of o — Ps — ;• invisible decay is the apparent disappearance of the energy 
lifig expected in ordinary decays in a hermetic calorimeter surrounding the o-Ps formation target. 
Therefore, the occurrence of the o — Ps^-o — Ps' — )• invisible conversion would appear as an excess 
of events with zero-energy deposition in the calorimeter above those expected either from Monte 
Carlo prediction of the background or from direct background measurements. 

The experiment presented here is based on the slow positron beam used to form o-Ps in a vac- 
uum cavity combined with the BGO calorimeter used in our previous search for o — Ps — )^ invisible 
decays [^l]], see also [ ]62| , 53 1. A preUminary design of the apparatus has been presented in Ref. 



54p. In the present work, more detailed simulations of the experimental setup, an improved 
positron tagging scheme, the selected positronium formation target and a better understanding of 
all possible background sources and sensitivity of the proposed search are presented. 

The great advantages of this approach compared to the one, where we produced o-Ps in an 



aerogel target [61], are listed here: 



Compared to our previous experiment, a factor 10^ more statistics can be collected with the 
same number of positrons. In the thin Si02 films that we plan to use as a target, 10 times 
more o-Ps is produced per implanted positron (see Section \\.2\ ) compared to the aerogel. 
Furthermore, there is no need to apply cuts for the 1.27 MeV photon selection and for the 
fiber energy deposition that reduced the number of events to less than 1% with respect to the 



number of positrons emitted from the source [61, 55]. 



• In the previous experiment, the main contribution to the 12% inefficiency for the detection 
of events that gave a trigger, arose from the overlap of close in time annihilation events (so 
called pileup events). In the beam based experiment described here, after every 300 ns bunch 
the chopper shuts off the positron flux for 3 /is. Therefore, the efficiency for signal detection 
will be close to 100% (see Section ^. 

• The suppression of the o — Ps ^ o — Ps — ;• invisible conversion due to the decoherence 
caused by the interaction of o-Ps with matter will be minimized. The number of collisions 
with the walls of the vacuum cavity undergone by the o-Ps during its lifetime is at least a 
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factor 10^ smaller than in the aerogel pores. Since the branching ratio for the o — Ps ^ 
o — Ps ^ invisible decay is approximately inversely proportional to the number of the o-Ps 
collisions (see Fig. ^, the sensitivity on the mixing strength will be a factor 100 better for 
the same statistics. 

In case an excess of events above the MC expectation for the background is observed, this 
experiment offers a unique and essential feature: one can cross check experimentally if it 
comes from signals. More precisely, by changing the number of collisions one can modify 
the oscillation probability while the background remains the same. We thought about two 
different possibilities: 



1. taking two runs at different positron implantation energies. From 3 to 5 keV the mean 
velocity of the created o-Ps increases by about a factor of two, thus, the collision rate 
with the walls is 2 times bigger and the signal is suppressed by the same factor (see 
Section O). 



2. the same result can be achieved varying the length of the cavity confining positronium 



(see Section |4.3[ ) while keeping the implantation energy of the positron fixed. 



However, compared to the previous experiment, there is a clear disadvantage: the calorimeter 
must be mounted outside the vacuum chamber so that the vacuum pipe introduces a loss of the 
photon energy. Nevertheless, simulations show that with an aluminum pipe of 1 mm thickness the 
sensitivity of the experiment will be at a level of 10^^ (see Sections ^ and ^. 



4. The setup 
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Figure 6. Schematic view oftlie setup. Tlie main parts oftlie system are: Tlie bundled slow positron beam, 
the secondary electron tagging system, the target chamber and the y detector 
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The experiment is designed with the goal to observe the o — Ps^ invisible decay if its branch- 
ing ratio is of the order of 10^^ (see Section ^. According to this requirement, the apparatus has 
several distinct and separated parts (Fig. H) that will be described in detail in the following sections: 



1 . the bunched slow positron beam (Section 4.1) 



2. the target for efficient o-Ps production near thermal energy (Section 4.2) 



3. the vacuum cavity to confine the o-Ps (Section |4.3| ) 



4. the positron appearance tagging system with a high S/N ratio, based on a high performance 
micro-channel-plate (MCP) described combined with the positron bunching (Section |4.4| ) 

5. the gamma detector, an almost An BGO crystal calorimeter (ECAL) surrounding the vac- 
uum cavity for efficient detection of annihilation photons to search for invisible o-Ps decays 



(Section 4.5) 



4.1 The slow positron beam 

The details and performance of the ETHZ positron beam are presented in [|6^. The beam is de- 
signed to operate in two modes, i.e. there are two different ways of tagging the positrons (or 
positronium): 

1. Detection with a MCP of the secondary electrons emitted when the positrons hit the target. 

2. Positron bunching: an initial pulse of 300 ns is compressed to the target region into a 2 ns 
wide pulse. 

In the first mode of operation, the secondary electrons (SE) produced by the positrons (about 25000 
e+/s) hitting the target are accelerated to 1-10 keV by the same voltage (applied to the target relative 
to the grounded transport tube) that is used to implant the positrons in the positronium converter 
(target). The secondary electrons are then transported in the backward direction by the same mag- 
netic field, which guides the positrons. The two operation modes described above can be combined. 
This is done by studying the delay between the detection of the SE at the MCP with respect to the 
arrival time of the positron at the target: the trigger for the positron tagging can be formed by a 
coincidence of the pulse from the MCP and the signal from the pulsed beam (see Fig. 0). This is a 
key point for the experiment because of the requirement to have the highest possible signal-to-noise 
(S/N) ratio in order to suppress fake triggers. With a fake trigger we mean that the trigger is not 
correlated with the presence of a positron in the o-Ps formation cavity. In this case, no annihilation 
photons would be detected and, as a consequence, this event would be mis-identified as an invisible 
decay. The requirement of this coincidence suppresses: 

• the background generated from electrons and ions produced by ionization of the residual gas 
atoms by positrons. 



the accidentals due to the MCP noise to a level of 10 



-5 
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This fundamental step of the experiment has been already tested successfully and the result is 
shown in Figure ^. As will be explained in Section 4^4 , to reach the required confidence level for 
the presence of a positron in the formation cavity, it is necessary to add an additional condition to 
the trigger scheme. 
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Figure 7. Timing between the secondary electrons detected in the MCP and the bunching pulse of the beam 
for 1 and 2 keV implantation energy of the positrons. 



4.2 o — Ps production target 



As reported in [57]-[^], systematic studies to characterize different samples that could be used 
as positronium converters for this experiment were performed. From those samples, we selected 



the so called F-samples (cf. [[70|]) for o-Ps production target of the experiment for the following 
reasons: 



1 . It provides a high production rate ofo — Ps in the range of 27-28 % for positron implantation 
energies between 3 and 5 keV (see Fig. HJ), thus, the high statistics required for the experiment 
can be reached (Section ^). 

2. In the interval of implantation energies between 3 and 5 keV the fraction of o-Ps remains 
constant within 1% (from 28% to 27%), thus, the background coming from the 2y decay 
oi p — Ps remains the same for the different implantation energies. While this fraction is 
almost constant, the mean energy of the produced o-Ps is about a factor of 3 larger for 3 keV 
implantation energy that for 5 keV (see Fig. ^. We estimated with the MC simulation that 
for a cylindrical cavity with 30 mm diameter and 15 mm length at 5 and 3 keV one has 0.9 
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and 1.9 collisions per lifetime. Therefore, the collision rate with the walls is ~ 2 times bigger 
for 3 keV positrons and the signal is suppressed by the same factor. 

It is worth noting that from 3 to 5 keV the implantation depth for the positrons varies from 
165 to 300 nm (the density of our porous films is about 1.5 g/cm^) thus the background will 
not be affected because the difference of the energy absorbed in 150 nm of porous SiOi is 
negligible. 

3. In the samples the o-Ps is produced near-thermal energy, thus, it experiences only few colli- 
sions per lifetime with the cavity walls minimizing the suppression of the signal. 

4. The target thickness is about 800 nm and it can be spin coated directly on the end plate of 
the beam pipe (Fig. ^) to avoid additional material of a substrate. Therefore, the absorption 
of photon energy in the target will be minimized. 
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Figure 8. (Top plot) Total fraction of o-Ps produced in the Si02 targets as a function of the positron implan- 
tation energy measured with positron annihilation lifetime spectroscopy. (Bottom plot) Emission energy of 
o-Ps in vacuum measured with time-of-flight (those data are taken from /p3/j- 



4.3 The vacuum cavity 

An essential issue in the experiment is that the positronium produced in the target should be con- 
fined in a region of highly uniform detection efficiency, and the o-Ps leakage through the aperture 
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where the positrons are implanted, has to be minimized. The o-Ps which escapes the detection 
region, could mimic an invisible decay because in this case the detection efficiency for the annihi- 
lation photons is reduced dramatically (Fig. ^ (a)). 
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Figure 9. (a) Ps is tagged by the secondary electrons produced by the positron that hits the converter 
surface (target). In the case where Ps escapes the detection region, it could mimic an invisible decay due to 
the strongly suppressed detection efficiency of the decay photons, (b) By closing the vacuum cavity around 
the o — Ps converter one obtains a region of highly uniform detection efficiency and the o-Ps leakage through 
the aperture is suppressed. The 15 nm carbon foil acts as a barrier for the o — Ps emitted from the converter 
but it is nearly transparent for the incoming positron beam and for the secondary electrons used for the 
trigger 



In their recent decay rate experiment ||75|], the Michigan group tried to minimize this effect to 
reduce the systematic error, using a double chamber-cavity, but still the escape probability for o-Ps 
was at 200 ppm. In our measurement such a leakage would limit the sensitivity to a Br(o — Ps — >• 
invisible) ~ 10^^. Therefore, to be able to reach the aimed sensitivity of the experiment, one has to 
invent a new method. We plan to completely close the o-Ps formation cavity ^ by employing a thin 
(15 nm) carbon film in which the few keV positrons can pass through. However, this will block 
even the most energetic positronium (several eV). Furthermore, with this method an additional 
signature for the presence of a positron in the formation cavity is added to the trigger scheme: the 
coincidence between the SE from the target and the ones produced in the carbon foil (Fig. ^ (b)). 
Hence, the S/N ratio will be further enhanced without a dramatic loss of the trigger rate. This 



point is presented in Section [4.4| , some preliminary measurements have been performed in order to 
test this idea. The carbon films that are planned to be used are currently employed at PSI"^ in the 



Only microscopic holes will be left for pumping the vacuum inside the cavity, leaving a negligible probability for 
the o-Ps to escape the detection region. 

"^PSI, Paul Scherrer Institut, CH-5232 Villingen-PSI, Switzerland 
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muonic -hydrogen experiment [ J77p for tagging the muons. Similar to our application they are used 
for detecting the SE emitted by the crossing muons. Some of those films were kindly given to us 
from the PSI group and we tested their SE emission yield in the implantation positron energy range 
of 1-7 keV and their permeability to o-Ps. The great advantage of this technique is that the o-Ps 
leakage is expected to be reduced to a negligible level. 

In order to minimize the energy loss of the annihilation photons in the vacuum cavity, one 
should design the beam pipe as thin as possible. In the HI experiment at DESY, a vacuum pipe with 
840 microns wall thickness was constructed. The internal part consisted of 40 microns aluminum, 
while the external part was made of carbon (800 microns). Similarly, we are designing the pipe for 
our experiment ^. 

To avoid problems with the magnetic field necessary to guide the positrons to the target region, 
the coil should be wound directly on the beam pipe, since we demonstrated with our TOE and PALS 



setup ([[78|]), that in this configuration the gain of the PMTs is not affected. Therefore, the material 
necessary for the coil should also be minimized. A coil made with a copper wire of 100 microns 
can produce the necessary field of 60 Gauss with a current of 0.5 A (the diameter of a wire fusing 
at this current is about 20 microns). In the following discussions, we consider a beam pipe with 
1 mm wall thickness, thus the expected sensitivity should be understood as a conservative result 
(Section |§). 

The carbon foil and its support could be suspended by some thin wires so that it is possible 
to move it in vacuum (see Eig. [l^. With such a scheme one could vary the distance between the 
target and the carbon foil, i.e., one would change the number of collisions that o-Ps will suffer 
during its lifetime. The signal will be suppressed for shorter distances in the same way as with the 
method described before by applying different positron implantation energies. 

4.4 Positron tagging system 



The positron tagging system (see Eig [11|) is based on the high performance MCP (Hamamatsu 
E4655-12) as a SE detector (noise rate < IHz). The coincidence between the SEs produced by 
positrons hitting the carbon foil and the SEs emerging from the o-Ps production target will be used 
to tag the positron appearance in the target region. 

The acceleration voltage for the positrons is applied in two steps. The positrons are first 
accelerated to 2 keV energy by the voltage applied between the ground pipe and the carbon foil, 
thus the SEs emerging from the carbon foil (CE) are transported with this energy back to the MCP. 
In the second acceleration stage, the positrons are accelerated by a 1-3 keV voltage applied between 
the CE and the target. Once they hit the target, the SEs are transported back to the carbon foil and 
the electrons able to cross it are then additionally accelerated so that their final energy will be 3-5 



keV (see Eigure [12|). Eigure |13| shows a MC simulation of the time distribution of the MCP signals 
from different SE sources; the time zero point is defined by the incoming beam positron crossing 
the carbon foil. The two prominent peaks are due to the SEs emitted when the positrons hit the 
target (at about 21 ns) and the SEs from the carbon foil (at about 29 ns) when a positron crosses 
the foil. The SEs from the target are accelerated with a higher voltage than those from the foil and 



^ We are also considering the possibility of making a vacuum pipe out of a scintillating material with the internal part 
coated with a thin aluminum foil. 
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Figure 10. Schematic of the target region and the carbon foil. 
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Figure 11. Schematic of the positron tagging system for the experiment using the MCP signal for secondary 
electrons emitted from the target and from the thin carbon foil. 



reach the MCP first, although they start later and have a larger flight path. The small third peak 
visible in the time spectrum of the MCP is produced by additional secondary electrons generated 
by some SEs from the target interacting in the carbon foil. Those will be accelerated by the same 
voltage that accelerates the positrons and, therefore, they are detected 1-1.5 ns after the second 
peak (the cavity length used in the simulation was 30 mm and the positron energy is 6 keV). The 
broadening of the first peak is due to the angular spread of the SEs emitted from the target and 
scattering in the carbon foil. 
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Figure 12. Trajectories of the positron (blue) and secondary electrons (red) for the new design with the 
carbon foil. 
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Figure 13. Simulated time distributions of the secondary electrons detected in the MCP. See text for details. 



4.4.1 Trigger efficiency and confidence level 

In order to estimate the trigger efficiency and fake trigger suppression of the tagging system de- 



scribed above, a test set-up has been constructed (see Fig. |14|). In this set-up the distance between 
the carbon foil and the target was 8 mm and the distance separating the carbon foil and the MCP 
was 90 cm. 

The measurements were performed by acquiring the MCP timing with a 5 Giga sampling 
oscilloscope (Lecroy wave-runner 44XMi). The threshold for the acquisition was set to 10 mV and 
a sample of 10^ events was recorded in a time window of 50 ns after the arrival time at the MCP 
of the first electron that was defined as START signal. The start signal can be produced by the SEs 
from the carbon foil or by the SE produced at the target. The probability to detect a SE coming 
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Ceramic isolating ring 




Figure 14. Setup used for the measurements with the tagging system. In order to accelerate the SEs from 
the target a voltage difference between the sample holder (green) and the carbon foil has to be applied. A 
ceramic ring electrically insulates the carbon foil holder and the target holder. 



from the target is suppressed by the angular spread introduced by crossing the carbon foil, when 
propagating in the backward direction. This introduces in the trajectories larger spirals along the 
magnetic field axis with respect to the SE from the CF, and thus, reduces the probability of the SE 
to hit the MCP active area. An example of an event recorded with the oscilloscope is shown in Fig. 
15. The START signal is given by SEs produced at the target and the pulse at 20 ns is from the SEs 



produced at the CF. The data in Fig. |16| show the time delay spectrum with respect to the START 
signal for the electrons arriving at the MCP with 10^ triggers. These results confirm qualitatively 
the prediction of the simulation^ . Note, that in this distribution the time is inverted with respect to 
the simulations. In the simulation the time t=0 is defined by the positron arrival time at the CF. 

The first peak at about 2ns (and a tail extending to 5 ns) is due to the MCP dead time: the first 
electron in the SE cloud from the target (or carbon foil) makes a START signal and after the dead 
time a second electron of the same cloud produces another signal. The second peak between 10 
and 17 ns is triggered by the SEs from the CF when the START signal is produced by SE from the 
target; only the events in this peak are used for the positron tagging. The broadening of this peak 
is due to the angular spread of the SEs emitted from the target and scattering in the carbon foil. 

The time separation of 15 ns between the START and the second peak is consistent with a 
simple estimation considering a straight propagation of the SE electrons produced at the carbon 
foil and the SEs from the sample (see Table IT]). 



Voltage [kV1 Flight time [ns] 



1.2 

3.8 



39 

22 



Af 17 ns 



Table 1. Estimated flight time for SE emerging from CF (1.2 kV) and from the sample 3.8 kV for a distance 
of 80 cm. The expected delay time is 17 ns. 



To estimate the trigger efficiency and the fake trigger suppression of our system, we compare 
two measurements: 



As will be discussed later in this section, Geant 4 does not reproduce correctly the transmission, scattering angles 
and backscattering coefficients of charged particles with few keV energies. 
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Figure 15. Example of an event acquired with an oscilloscope. The first signal is produced by the SE from 
the target and defines the START signal. The second pulse is produced by the SEfrom the CF ^ I5ns after 
the START signal. The time window O) between 15 and 22 ns (i.e. between 10 and 17 ns delay from the 
START signal) is defined in order to calculate the trigger efficiency and the fake trigger ratio. A Gaussian 
fit is used in order to determine the position of the pulses (green). The origin of the third peak are the SE 
produced by the SEfrom the target crossing the CF as expected from the simulation (see Fig. 13). 
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Figure 16. Measured delay of the SEfrom the carbon foil relative to SEfrom the START signal. The dashed 
line represents the data obtained by applying the same voltage on the carbon foil and on the sample (-1.2kV). 
The solid line is obtained by applying -1.2 kV on the carbon foil and -3.8kV on the sample. In this case the 
broadening of the peak at around 15 ns is mainly due to variations of the fiight time of the SEfrom the target 
that produce the START signal. The peak at about 2 ns are due to the MCP deadtime (see text for details). 
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a set of data with a voltage difference between the carbon foil and the sample of -1.2 kV and 



-3.8 kV, respectively (solid line in Fig. 16), 



a set of data with the same voltage applied on the carbon foil and the sample (-1.2 kV) 



(dashed line in Fig. 16). 



We define the trigger efficiency Ese as: 

^SE = ^ (4.1) 

where Nw is the number of events arriving with a time delay from the start signal in a selected time 
window w when the a voltage difference is applied between the target and the carbon foil. Ntot is 
the total number of START signals. 

The carbon foil suppression factor Kse is defined as: 

n' 

Kse = -^ (4.2) 

where A^,^, is the number of events in the selected time window w when the same voltage is applied 
to the sample and the CF (i.e. SE produced at the target cannot reach the MCP). In this case the 
probability for SEs emitted from the sample to reach the MCP is strongly suppressed because the 
transparency of the foil for those electrons is practically zero. The START signal and the peak 



at ~ 2 ns for the dashed line in Fig. |16| (the same voltage on target and CF) is produced by SEs 
from the carbon foil, emitted by an incoming positron, since (almost) no SE from the target reach 
the MCP. Thus, Kse gives suppression of fake triggers obtained by the coincidence positron trigger 
requirement for events with a positron traversing the carbon foil, but no SE emitted from the target. 
It is used to estimate the background from positrons, which do not reach the target, e.g. due to 



backscattering from the carbon foil. This background will be discussed in Sections 5.1 and 5.2. 
Note, that the value Kse (see below) is not in contradiction with the requirement for the invisible 
decay search experiment to have a fake trigger fraction with no energy deposition in the ECAL of 
< 10^^. The suppression Kse refers to events with a positron passing the carbon foil and for most 
of these events the annihilation photons will deposit energy in the ECAL. 

The signals of SEs emerging from the sample and those emerging from the CF can be further 
discriminated by their charge. The large scattering angle acquired by the SEs from the target 
crossing the CF drastically reduces the probability to have, as a trigger, a signal produced by more 
than one electron. On the other hand, the STOP signal is composed of a larger number of SE 
because it is produced by the sum of: 

• the SE emitted by positrons crossing the CF; 

• the SE that are emitted when, the SE produced at the target, cross the carbon foil in the 
opposite direction, since the small difference in their arrival time at the MCP (relative to 
the SE from the oncoming positron) is not resolved because of the small distance (8 mm) 
between the target and the carbon foil and the large voltage between them. 
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This is clearly visible in the Fig. |17| where the charge of the START signal (i.e. the first 
electrons which hits the MCP) is plotted versus the signal in the time window of 12-17 ns for both 
voltage configurations. This effect can be exploited to considerably increase the confidence level 
for a positron in the vacuum cavity without dramatically affecting the trigger efficiency as shown 



in Fig. 18. The trigger efficiency (boxes) and the coincidence suppression factor k:^^ (circles) are 



studied for two different selection cuts on the MCP signal: 

1. by applying a threshold on the MCP stop signal; 

2. by applying a cut on the charge (integral over the signal) produced by the MCP. 

In the first case the fake trigger ratio cannot be suppressed below a value of 4x10^^ while the 
trigger efficiency drops to a value of 10^^. In the second case, the cut performed on the charge 
shows to be much more efficient. For a cut on the charge of 40 pVs is suppressed to 5x 10^^ with 
a corresponding trigger efficiency of 4x 10^^. 

4.5 The photon detector 

The same BGO crystals as in our former o-Ps invisible decay search will be used to detect /-quanta 
produced in positron or positronium annihilation. The geometry of the detector should be modified 



to accommodate the beam pipe as we proposed in [64] (see Fig. 19). The test of the trigger 



scheme presented in Section i.4 was done in a straight geometry, however the end of the beam pipe 
should be bent in order to avoid a direct way for the annihilation's photons and for Ps to escape 
the detection region. We are not expecting a decrease of the trigger efficiency or an increase in the 
background because our simulation predicts that with a careful design the transportation efficiency 
for both positrons and SE through the curved beam pipe region is 100%. Furthermore, the crystals 
have to be arranged perpendicularly to the B -field, since, as we confirmed in our TOF and PALS 



setup ([|78|]), in this configuration the magnetic field will not affect the performance of the PMTs. 

The 7-detector serves to veto effectively the positron annihilations into photons. It has been 
shown that its inefficiency for the detection of annihilation's photons is less then 10^*^ for the en- 
ergy threshold of ~ 80 keV ^]\ , ^^. Moreover, we realized that the energy resolution of the 
crystals could be increased by replacing the Teflon in which the crystals are wrapped with the 3M 
radiant foil. This reduces the amount of dead material by a factor 12. Furthermore, all the BGO 
crystals surface is now polished. By roughing the BGO surface except where the PMT is coupled 
to the crystal, the energy resolution of the calorimeter could be further improved. Therefore, the 
calorimeter could be refurbished in order to increase the energy resolution resulting in an improve- 
ment of the detector performance. 

5. Background estimation 

As mentioned in the previous section, we expect backgrounds which originate from the following 
sources: 



^It is worth noting that in the experiment described in this section the threshold will be lower because the pileup will 
be reduced by a factor 20. 
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Figure 17. Charge (measured with 50 D, termination) in pVs of the start signal versus the charge of the 
signal for events in the time window 12-17 ns. Top: the same voltage of 1.2 keV is applied to the target and 
the carbon foil. Bottom: a voltage of 3.8 keV is applied to the target and 1.2 keV to the carbon foil. 



1. the annihilation energy losses that are estimated to be at a level of ~ 10^^. 

2. a fake positron tagging at a level ~ 10^' 

The main contribution to the first background of the list above is coming from the thickness 



of the vacuum beam pipe. In Fig. 20, we show the simulation results of the energy deposited in 
a 0.84 mm thick aluminum pipe and a pipe with 0.04 mm aluminum and 0.800 mm thick carbon 
pipe (similar to the one that was used at the HI experiment at DESY). In these distributions, the 
energy deposited in the target substrate and in the copper wire surrounding the beam pipe are also 
included. 

Different possibilities of fake triggers that can be produced in this experiment have been iden- 
tified: 
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Figure 18. (Top plot) Trigger efficiency (boxes) and coincidence suppression factor KgE (circles) as a func- 
tion of the threshold set on the signal in the time window 12-17 ns. (Bottom plot) Trigger efficiency (boxes) 
and coincidence suppression factor K^e (circles) as a function of the cut on the charge (voltage integral) on 
the signals in the time window 12-17 ns. 



• Fast backscattered o — Ps produced at the carbon foil surface have a certain probability to 
escape the detection region. 

• Backscattered positrons from the carbon foil surface or the target have a probability to escape 
the detection region. 

5.1 Fast backscattered o — Ps from the carbon foil 

Fast backscattered o — Ps can be produced from shallow implanted positrons that capture an elec- 
tron exiting the surface. The contribution to the background arising from this effect has been 
studied with the MC simulation. We approximated the energy distribution of o — Ps with a Landau 
distribution peaking around 15 eV [ [7T[ ] (see Fig. ^1|)- The probability P{E) of o — Ps to annihilate 
via pick-off when it collides with the cavity walls was also included. As a function of the energy, 



P(E) can be divided in two regions ([76]): 



1. For Eo-Ps < 6.?,eV , P{E) = 

2. For Eo^p, > 6.SeV , P{E) = 0.95 
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Beam pipe 




Figure 19. Cross section of the BGO calorimeter mounted around the beam pipe. The sphere is to show that 
the minimal BGO thickness around the target is of the order of 200 mm. 



Both assumptions are conservative. In fact, even o — Ps with an energy smaller than its binding 
energy (6.8 eV) can undergo pick-off annihilation, thus the number of Ps atoms escaping the detec- 
tion region will be less than what we estimated with the simulation using the assumption 1). The 
same is true for Ps with energies above its bounding energy because in this case the probability to 



dissociate in a collision is close to 100% [76]. We define the escaping probability t, as: 



^o-Ps 



N, 



escape 
Nrnr 



(5.1) 



where Nescape is the number ofo — Ps which decay after the bending of the vacuum cavity (See Fig. 
22) and Ntot is the total number of events simulated. Note that this approach is very conservative 
because the annihilation gammas produced outside the bending region have a non-zero probability 
to deposit some energy in the calorimeter. The escaping probability estimated with the simulation 
(Fig. 22) is t,o-ps ci 1 X 10^"*. In order to calculate how this value affects the sensitivity of the 
experiment, it has to be multiplied with the fast o — Ps formation probability [[7T[] (< 10%). In 



order to mimic an invisible decay, the fast o — Ps escaping event has to coincide with an accidental 



trigger. Thus, the carbon foil suppression factor Kse presented in Section [4.4. 1| further suppresses 
this background to a level < 5 x 10^^. 

5.2 Backscattered positrons 

The positrons are transported along the beam line with an energy that can be varied from Eq = 
10 — 200 eV. The carbon foil and the target are biased with a potential with respect to the beam 
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Figure 20. Distributions of the energy deposited in the dead material surrounding the target region (target 
substrate, beam pipe and copper coil) from annihilation events in the target. The top plot shows the results 
of the MC simulation for a 0.84 mm thick aluminum pipe. The bottom plot the distribution for the pipe 
construction with 0.04 mm aluminum and 0.800 mm carbon. The total number of simulated 2y-events is 10 
in both cases. The peaks at 511 keV and 1022 keV correspond to the total photo-absorption either of a single 
511 keV photon or of both of them, respectively. 
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Figure 21. Energy distribution of fast o-Ps originating from backscattered e^. The pick-off probability after 
a collision of o-Ps with the walls is taken to be zero for o-Ps energies smaller than the 6.8 eV while for 
higher energies it is 0.95. 



pipe, giving an additional acceleration to the positron, thus, their energy at the carbon foil is Etot = 
Eq + Ecf and at the target it is Ejot = Eq + Ejarget- If a positron backscatters in the carbon foil 
(or in the target) it is decelerated by the applied voltage. In this case, it can either be re-implanted 
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Figure 22. Simulated decay position ofo-Ps in the vacuum cavity. The escaping probability estimated with 
the MC simulation w ~ 1 x 10^^. 



or escape the detection region. The fate of the backscattered positron depends on its energy loss in 
the backscattering process at the foil or target AE = Ejot — Esack'- 

1 . if EBack < EcF (Target) the positron is re-implanted in the carbon foil (or in the target) with an 
energy EBack- 

2. if Esack > Ecpijarget) the positron has enough energy to escape the carbon foil (the target) 
voltage and it can be transported back in the beam line with an energy Esack — EcF(Target) 

This second process represents a dangerous source of background. An escaping backscattered 
positron could give a trigger and be transported back by the magnetic field outside the detection 
region, and thus, it will not deposit any energy in the calorimeter. This effect can be suppressed by 
decreasing the beam transport energy £"0 so that the minimum energy loss necessary to escape from 
the target voltage is smaller and consequently also the escaping probability is suppressed. In this 
section, a precise estimation of this effect will be carried out using MC calculations. The dominant 
process in the backscattering effect is the positron multiple-scattering. GEANT4 has a set of low 
energy classes that can be included in the Physics list using PENELOPE MC cross sections such as 
G4PenelopeBrehmsstrahlung, G4PenelopeComptonEjfect. Unfortunately, the multiple-scattering 
process is not yet implemented and consequently it is not possible to estimate the backscatter- 
ing coefficients with GEANT4. Accurate knowledge of the energy and angular distributions of 
keV charged particles is possible with the EGSnrc (Electron Gamma Shower) Monte Carlo code 



[72]. This code includes sophisticated low-energy physics comparable with PENELOPE. More- 
over, EGSnrc is considerably faster than PENELPOPE and does not require careful tuning of the 
simulation parameters. In EGSnrc a special package was developed for the calculation of backscat- 
tering coefficients. The performance of EGSnrc for keV backscattering particles is presented in 



[72] where MC simulation data are compared to experimental results for both electron and positron 



backscattering (see Fig. |23| ). 
We define the positron escaping ratio Eg+^esc from the detection region as: 
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Figure 23. Positron backscattering coefficients (ri+) versus positron kinetic energy Iy2\ 



■-NeJNu, 



(5.2) 



where Ne„ is the number of backscattered positrons whose longitudinal energy loss is smaller than 
the incoming positron energy before the acceleration. 

The effect of the positron backscattering has been studied in two situations: 

1. Positron backscattering at the carbon foil. 

2. Positron backscattering at the SiOi target. 

5.2.1 Positron backscattering from the carbon foil 

To calculate the backscattering probability, we assume that the carbon foil has a density of Igjcnv' 
and a thickness of 20 nm. There is a substantial difference between the backscattering coefficients 
for semi-infinite targets and thin foils because of the transparency effect that shows a considerable 
suppression of the backscattering probability when the incident energy of the positron increases. 
In Fig. ^ the backscattering coefficient ?]+ for a thick carbon target and for a 20 nm carbon foil 
are compared. On the same plot the transparency of the foil for the same implantation voltages is 
shown. 
The probability of a positron to escape the carbon foil voltage depends on: 

• the backscattering coefficient 77+ 

• the initial positron energy E^ before the foil acceleration Eqf 
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the backscattering energy and angular distributions 
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Figure 24. Positron backscatter coefficient (rj+jfrom a thicJc carbon target (black squares) and a thin carbon 
foil (triangles) as a function of the positron kinetic energy. The transmission probability, through the thin 
foil is also shown (open squares, right scale). 



In order to estimate the escaping fraction, the whole energy distribution of the backscattered 



positrons has to be known (see Fig. 25). Figure 26 shows the escaping probability for positrons 
backscattered from the carbon foil as a function of the acceleration voltage. The calculations have 
been performed for different positron initial energies Eq. As one can see, for a good choice of the 
voltage applied to the carbon foil (about 1.5 kV) and of the initial energy (the lower the better, 
Eq = 10 eV) the escaping probability is smaller than 2x10^^. Positrons backscattering at the 
carbon foil are not a direct source of background. To result in an artificial invisible decay, such 
an event has to coincide with an accidental trigger. Therefore, this background is suppressed by 



the confidence level of the tagging system presented in Section 4.4.1 and results in a background 
< 10^. 



5.2.2 Positron backscattering from the Si02 target 

After crossing the carbon foil the positron enters the cavity and it is accelerated to the Si02 target 
where it may backscatter. For this estimation, we consider the S/O2 target density of Ig/cm^. The 
backscattering coefficient for Si02 is larger than the one of the 20 nm foil (Fig. p^. In agreement 



with the results reported in [73|, 77+ reaches a constant value of about 10% for positron incident 
energies larger than 4 keV. The positron escaping probability as a function of the energy is shown 
in Fig. ESl For positron energies between 3-5 keV the escaping probability is between 1 x 10^^ and 



3x 10^^. These values are larger than the maximum background level that is allowed in order to 
reach the desired sensitivity. However, this background is clearly overestimated because the energy 
loss of the positron crossing the carbon foil is not taken into account. An estimation of the mean 



energy loss [74] [72] gives values of the order of few hundred eV that are significantly larger than 
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Figure 25. (Top plot) simulation using EGSnrc. Energy spectrum of positrons backscatte red from a 20 nm 
carbon foil for positron implantation energies in the 1-2 keV range. The cut-off energy is 300 eV. (Bottom 
plot) angular distribution for backscatte red positrons 
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Figure 26. Escaping probability for positrons backscatte red from the carbon foil as a function of the carbon 
foil voltage for different positron initial energies Eq. 



the initial positron energy £0 = 10 eV. 

A possible way to further suppress this background is to redirect to the target the positrons 
that may backscatter. The idea is to exploit the difference in the propagation time between the 
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Figure 27. Backscattering coefficient r\+for the SiOi target as a function of the positron incident energy. 
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Figure 28. Escaping probability for the total number of incident positrons at the Si02 target as a function 
of the implantation voltage for different positron initials energies Eq. 



backscattered positrons (< 10 eV) and the SE used for the trigger (>2 keV). When a SE triggers the 
MCP an electrode with a repulsive voltage can be activated to block and redirect back a positron 



that may have backscattered as shown in Fig. g9[ The optimal position of the repulsion electrode 
has been determined assuming that the backscattered positrons have an energy between 0. 1 and 
10 eV. We defined the minimum distance between the electrode and the target considering the 
propagation time of the SE from the target to the MCP (70 ns with d=l m, E= 2 keV) and the 
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propagation time of the backscattered positron. The maximum electrode-target separation is given 
by the constraint that the detection probability for a Ps atom that is eventually formed when the 
positron is re-implanted in the target should still be high enough to guarantee the aimed sensitivity 
of the experiment. In fact, the propagation time of the positron from the target to the repulsion 
electrode and back shorten the acquisition time window for such an event. Taking into account the 
backscattering probability and the Ps lifetime, we estimated that one needs an acquisition time of 
at least 0.5 /is. With these requirements the electrode should be positioned in the region indicated 



in Fig. 29 with a thick line (green). 
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Figure 29. Scheme of the repulsion electrode used to suppress the backscattered positrons and plot with the 
optimal position (thick green line on the y-axis). See text for more details. 



In this way, the background produced by backscattered positrons from both the carbon foil or 
the target can be suppressed to values below < 10^'. Table^ summarizes the expected background 
level for the different background sources. 

6. Sensitivity 

The sensitivity ^o-ps^invisibie of the experiment is defined as the level at which the first background 
event is expected: 

5'o-Ps^invisible = '^/{No-Ps ' £tot) (6. 1) 

where the terms in the denominator are the integrated number of produced o — Ps (No-Ps) and Etot — 
1 is the total efficiency to detect an invisible decay. We neglected the losses in signal efficiency 
of about 1.5% arising from the possibility of having 2 or more positrons per bunch. We estimated 
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Table 2. Summary of the expected background level for the different background sources. 



the rate of these events using /?2e+ = 2 • Thunch Re+ where Tbunch = 300 ns and Rg+ = 2.5 x W^ /s is 
the number of delivered positrons per second on the target in continuous mode. For two or more 
positrons there is always annihilation energy deposition in the ECAL, hence this effect does not 
result in a background. The number of o-Ps/s, Ro-Ps, is defined as a product 



J^o—Ps — ^e+ ■ ^o—Ps ' ^SEtagging ' ^SEcoincidence ' ^Bunching 



(6.2) 



where the first factor was defined above, the second one is the efficiency for o — Ps production 
(about 30%) and the third one EsEtagging is the efficiency of tagging the secondary electrons from 
the carbon foil (about 20%). The fourth factor is the efficiency of the coincidence between the 
SEs from the CF and target {esEcoinddence = 0.04) and the last one, SBunchmg = 0.1, are the losses 
due to the duty cycle of the bunching system. As in our previous search, the length of the gate 
for the ADCs has to be at least 3 /is in order to suppress the probability for o-Ps to decay after 
this time to a level of 10^^. Therefore, the aimed sensitivity of 10^^ can be reached in a 8 hours 
run (rs 1 X 10^ observed o — Ps annihilations). For zero signal events observed and no event of 
background expected, the upper limit at 90% CL for the branching ratio assuming Poisson statistics 
is given by: 

Br{o -Ps^ invisible) = 23/{No-Ps ■ £tot) (6-3) 



Solving Eq. (2. 17 ) with the estimated average number ofo — Ps collisions in the cavity Ncoii — 0.9 
(for 5 keV implantation energy of the positrons), shows that a limit on the Br{o — Ps — ;• invisible) = 
1 X 10^^ will result in a limit on the photon-mirror photon mixing strength of: 



£<4xl0"^ 



(6.4) 



This is about one order of magnitude more stringent than the BBN limit ( |2.21[ ) and can be achieved 
in a 1 8 hours run. 
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Assuming that the DAMA/Nal and DAMA/LIBRA annual signal modulation is generated by 
elastic scattering of mirror matter, the mixing strength is of the order of e < 4 x 10^^, thus a total 
number of ~ 93 signal events would be detected in the ECAL during one month of data taking. 
To be conservative we consider a background level of 1 x 10^^, thus, about 90 background events 
are expected which means that a discovery with about 7 a significance could be possible [p4|]. As 
explained in Section ^ a unique feature of our proposal is the possibility to change the experimental 
conditions (i.e. the number of the o-Ps collisions with matter), and hence to cross check the results 
without affecting the background. For an implantation energy of the positron of 3 keV, the number 
of excess events will be 2 times smaller (46 events) compared to 5 keV positrons. 

7. Summary 

In this paper, a proposed search for mirror-type dark matter, looking for o-Ps invisible decays in 
a vacuum cavity, is presented. In the Introduction, the mirror matter relevance to the dark matter 
problem and its link with positronium has been reviewed. The effect of external fields (electric and 
magnetic) on the oscillation probability and of the collisions of o-Ps with the cavity walls were 
estimated. In Section ^ the design and the experimental results for the different parts of the exper- 
iment have been presented. Section ^ includes the simulation results of the background estimation 
and Section ^ the estimated sensitivity. The goal is to reach a sensitivity in the branching ratio of 
Br{o — Ps — )• invisible) ~ 10^^ to confront the annual modulation signal observed by DAMA/Nal 
and DAMA/LIBRA (with 8.2a significance) with mirror dark matter scenarios. In case of a signal 
observation, the experiment would offer a unique and essential feature, whereby an increase or 
decrease of the signal rate by a factor ~ 2 is possible while keeping the background constant. In 
case of a null result, this search will provide a limit on the photon-mirror photon mixing strength 
about one order of magnitude better than presently derived from Big Bang Nucleosynthesis. 
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